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ABSTRACT 

We examine scatter and bias in weak lensing selected clusters, employing both an 
analytic model of dark matter haloes and numerical mock data of weak lensing cluster 
surveys. We pay special attention to effects of the diversity of dark matter distributions 
within clusters. We find that peak heights of the lensing convergence map correlates 
rather poorly with the virial mass of haloes. The correlation is tighter for the spherical 
overdensity mass with a higher mean interior density (e.g., Miooo)- We examine the 
dependence of the halo shape on the peak heights, and find that the root-mean-square 
scatter caused by the halo diversity scales linearly with the peak heights with the 
proportionality factor of 0.1 — 0.2. The noise originated from the halo shape is found 
to be comparable to the source galaxy shape noise and the cosmic shear noise. We find 
the significant halo orientation bias, i.e., weak lensing selected clusters on average have 
their major axes aligned with the line-of-sight direction, and that the orientation bias 
is stronger for higher signal-to-noise ratio (S/N) peaks. We compute the orientation 
bias using an analytic triaxial halo model and obtain results quite consistent with 
the ray-tracing results. We develop a prescription to analytically compute the number 
count of weak lensing peaks taking into account all the main sources of scatters in peak 
heights. We find that the improved analytic predictions agree well with the simulation 
results for high S/N peaks of > 5. We also compare the expected number count with 
our weak lensing analysis results for 4 deg^ of Subaru/Suprime-Cam observations and 
find a good agreement. 
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INTRODUCTION 

' Clusters of galaxies have been playing an important role 
in the field of cosmology. For instance, number counts of 
clusters have pl aced useful constraints on cosmological pa- 
rameters (e.g., Vikhlinin et~all l2009bl : iMantz et all I2OI0I : 
IRozo et al.ll201C )■ and measurement of the dark matter dis- 
tributio n in clusters is useful for testing d ark mat ter mod- 
els (e .g., 
20081 : 
20121: 



cosmic microwave ba ckground (e.g., [Marriage et all I2OIII : 
iReichardt et all I2012I '). and dark matter concentrations 
in the weak lensing mass ma p (e.g., ^Miyazaki et al 



,2002, 2007; W ittman ct al 2006,; iGavazzi fc Soucail 2007 
Schirmer et al.! l2007l: [Kubo et al.l l2009l: iBeUagamba et all 



Dahlell2006l:lBrqadhurst et al.|[2008 : Mahd avi et"all 



200C 



Okabe et al. 



20101. 



Okabe fc Um etsu' '2008*: 'Oguri ct al. 
Umetsu et al . 2009, 2011; Mcdczinski ct aj] 

2OIOI) . In most of these apphcations, the con- 



20ld : 



struction of homogeneous cluster samples and understand- 
ing selection biases are of fundamental importance, because 
an unknown bias can significantly affect the interpretation 
of the data. 

Clusters of gala:xies are identified by various techniques, 
including detections of galaxy con centrations in optical 
data, extended X-ray emissions fe.g.. lBohringer et al. I I2OO4I : 
I Vikhlinin et al.ll2009al ). the Sunyaev-Zel'dovich effect in the 



2OIII : IShan et al.l I2OI2I : iKurtz et all 120121 1? Each method- 
ology has its own advantage and disadvantage. Among 
others, weak lensing technique is unique in the sense 
that it identifies clusters by searching for high peaks in 
the weak lensing mass map, and thus d oes not rely on 
physical state of the baryoni c compo ne nt (ISchneideilll996l: 
Harna,na. Takada. fc Yoshidal |2004| : iHennawi fc Spergdl 



20051 : iMaturi et al."2005''). Therefore, combining weak lens- 



ing cluster sample with samples selected by another method 
is beneficial not only for various cosmological applications 
but also for better understanding of cluster physics. 

Recently IShan et al.l (|2012l ') reported 301 weak lens- 
ing high peaks located from 64 deg^ data. Among those 
peaks, they confirmed 85 groups/clusters, which is the 
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lars est weak lensine cluste r sa mples constructed to d ate (see 
also lMivazaki et al.ll|2007h and lSchirmer et al.l (|2007l )'). Cur- 
rently the size of weak lensing selected cluster sample is 
mainly limited by areas of optical imaging surveys. How- 
ever, the situation is changing drastically in coming decade 
because there are several ongoing or planned wide field deep 
optical surveys including Dark Energy Surve^Q, the Kilo- 
Degree Surveys, Hyper Suprime-Cairlj, and Large Synoptic 
Survey Telescop4f|- Thus in the near future, weak lensing 
selected cluster catalogues containing ~ 10, 000 clusters will 
be available. It is therefore important to examine the selec- 
tion function of the weak lensing selected cluster samples 
in order to take full advantage of these unique cluster cata- 
logues. 

The standard method for identifying clusters with weak 
lensing is to search for high peaks in the weak lensing mass 
map generated from weak lensing shear data with a carefully 
designed smoothing filter. Peaks above a given threshold are 
selected as cluster candidates. The peak height, which plays 
a central role in this study, is primarily determined by the 
mass and rcdshift of the da rk matter halo of clusters (see 
e.g. , Hamana, T akada. fc Yo shida 200i) , but is also affected 
by the following three effects. The first is the noise arising 
from intrinsic shapes of source galaxies used for weak lensing 
shear measurements (gala:xy shape noise). The second is the 
projection of structures along the line-of-sight. The third 
is the diversity of dark matter distributions in individual 
haloes. These effects are important in the sense that they 
induce the scatter and bias in the peak heights and hence 
in the resulting weak lensing selected cluster sample. While 
these effects were examined in literature to some extent 
(IHamana. Takada. fc Yoshida [2004: Hennawi fc Spergell 



flMamana. iakada. fe Yoshidal 12004: Mennawi fe spergel 
20051: iMaturi et al II2005I: iTan fc FanlfeoOSi: iPace et al.l | 2007 



Fan, Shan fc Liu"20ld : iMarian. Smith, fc Bernsteinl boiC ; 
Schmidt & Rozo 2011: iDietrich et alj|2012h . we revisit this 
problem with a particular emphasis on the scatter and bias 
produced by the diversity of dark matter distributions, 
using both analytic and numerical approaches. 

The matter content of clusters is dominated by dark 
matter, which accounts for ~ 80 per cent of the mass of 
the universe. The dark matter distribution within dark mat- 
ter haloes are investigated in detail using A'^-body simula- 
tions, which found that the radial profile of the dark matter 
halo is accurately describ ed by an a nalytic form, the so- 
called NFW profile (,Navarro. Frenk. fc White 1997 '). In this 
model, the matter distribution of haloes is characterized by 
two parameters, the mass and the concentration parameter. 
While it is known that there is a mean relationship between 
the mass and concentration parameter, the re lation involves 
a lar ge scatter among diff erent haloes (e . g.. iBuUock et al.l 
120011 ). As pointed out by iKing fc MeadI l|201ll ). the weak 
lensing peak height is sensitive to the halo concentration. 
Another important prediction of the current standard struc- 
ture formation model is that the dark matter distribution 
in clus ters is not spherica lly symmetric but is highly elon- 
gated (jjing fc Sutdl2002h . This non-sphericity of the dark 



^ http : //www . darkenergysurvey . org/ 
^ http://www.astro-wisc.org/projccts/KIDS/ 
^ http: / /subarutelescope . org/Projects/HSC/ 
* http://www.lsst.org/lsst/ 



matter distribution is known to h ave large impact on weak 
lensing measurements of clusters ( Clowe. De Lucia, fc King 
'2004'; Qguri ct al.' 2005'; 'Cavazzi 2005': 'Tan fc FanI l2005l : 
Corl css fc King. 2009 : . Feroz fc Hobson 2012). Thus even the 
same mass haloes at the same redshift can produce largely 
different peak heights depending on the halo concentration 
and non-sphericity. This is exactly what we explore in this 
paper using a large set of mock data generated from gravi- 
tational lensing ray-tracing in TV-body simulations. 

The structure of this paper is as follows. In Section[2]we 
summarize the analytic description of cluster density distri- 
bution adopted in this paper and also summarize basics of 
cluster finding with weak lensing. We use a mock numerical 
simulation of weak lensing cluster survey which is detailed in 
Section[3l In Section|4]we investigate the scatter and bias in 
the weak lensing selected clusters, paying special attention 
to the infiuence of halo shape. We briefly compare the peak 
number count derived from the numerical mock data with 
our weak lensing analysis results of Subaru/Suprime-Cam 
observations in Section [5] Finally, summary and discussion 
are given in Section [6l 

Throughout this paper we adopt the cosmological model 
with the matter density Q.m = 0.238, baryon density fij, = 
0.042, cosmological constant 0,^ ~ 0.762, spectral index 
Tig = 0.958, the normalization of the matter fluctuation 
(Ts = 0.76, and the Hubble parameter h = 0.732, which 
are the best-flt cosmological parameters in the Wilkinson 
Microwave Anisotro py Probe (WMAP) third-year results 
ISpergel et al.ll2007l '). 



2 ANALYTIC MODELS OF DARK MATTER 
HALOES 

Let us flrst summarize the analytic models of dark matter 
distribution of clusters of galaxies and its weak lensing prop- 
erties. 



2.1 Dark matter distribution within haloes 

We adopt the NF W model for the dark halo dens ity proflle 
which is given by (jNavarro. Frenk. fc Wh"it3ll997f l. 



PNFw(r) 



(r/rs)(l + r/r,)2- 



(1) 



This model characterizes the dark matter distribution by 
two parameters: the density parameter ps and the scale ra- 
dius rs- It is customary to re-characterizes it by related pa- 
rameters, the halo mass and the concentration parameter, 
which we deflne below. 

The mass of haloes is not a uniquely defined quantity. 
One needs to choose proper definitions of masses depending 
on the purposes (White 2001). In theoretical and observa- 
tions studies, the virial mass Mvir, which is defined such 
that the average density within the virial radius becomes 
equal to the nonlinear overdensity Ayjr computed using the 
spherical collapse model (see, e.g.. iNakamura fc Sutolll997t ) 
times the mean matter density of the universe, has often 
been adopted. For the NFW model, the virial mass relates 
to ps and Ts by 

47r 2 3 

Mvir = — Avir(z)pm(z)rvir = 47rp3r^ rTlnfw (Cvir ) , (2) 
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where Cvir is the so-caUed concentration parameter defined 
by 

1/3 

(3) 



1 



47rAvir(2:)pm(2:) 

and mnfw(cvir) defined by 

X 



mnfw(Cvir) 



(1 + X)2 



dx = ln(l + Cvir) — 



Cvir 



1 + Cvir 



(4) 



The concentration parameter is known to be correlated with 
the halo mass and redshift. When necessary, we adopt the 
following relationship: 

-0.086 



Cvir(Mvir,2) = 7.26 



Mvir 



1012/i-iMq 



(1 + ^)" 



(5) 



which was derived from A'^-body simula- 
tions assuming the WMAP third year results 
iMaccio. Dutton. fc van den BoschI I2OO8 ). with the ad- 
di tional redshi f t dep endence based on the simulation result 
of iDuffv et all ||2008D . 

While the virial mass is a physically motivated defini- 
tion of the halo mass, one can define halo masses using ar- 
bitrary values of overdensities. Specifically, one can use the 
spherical overdensity mass Ma defined by the mass con- 
tained within a radius ta inside of which the mean interior 
density is A times the critical densitjjf] 

Ma = — Apcr(z)rA = 47rpsr^mnfw(cA), (6) 
where 



CA 



rA 



3Ma 



-,1/3 



47rApcr(z) 



(7) 



In this paper, w-e also consider a triaxial halo model 
of lJing fc Sutol (|2002l ') to investigate the effect of the halo 
triaxiality. In this model, the density profile given by eq. ([1]) 
is modified as 

Pcd 



Ptri(r) = 



7?" 



(R/Ro){l + R/Ro) 



2 2 2 

X V Z , , , . 



(8) 



(9) 



IJing fc Sutol (|2002h derived the probability distribution of 
axis ratios a/c and b/c for haloes with a given mass and red- 
shift from large A'^-body simulations. T he lensing properties 
of the triax ial halo model were deri ved in lOguri. Lee, fc Sutol 
l|2003t ) and lOguri fc KeetonI (|2004l '). including projections of 
triaxial haloes along arbitrary directions, which we adopt in 
this paper. We study the effect of the halo triaxility on the 
properties of weak lensing selec ted clusters using the sem i- 
analytic approach developed bv lOguri fc BlandfordI l|2009l ). 
In this method, we generate a catalogue of haloes according 
to the mass functi on as well as axis ratio distributions of 
jjing fc Sutol (|2002i ). and project each halo along random di- 
rection to compute the lensing properties. This allows us to 
generate a mock catalogue of weak lensing selected clusters 
based on the triaxial halo model. 

^ Some authors adopt a different definition in which the A is 
specified relative to the average background density prn{z). Our 
A is Qrn{z) times theirs. 



2.2 Basics of weak lensing cluster finding 

Here we summarize equations for weak lensing clus- 
ter finding which are directly relevant to the fol- 
lowing analyses. For more detailed descriptions, see 
Schneider (1996); Bartelmarm, King_fc Schneider (2001 



Hamana. Takada. fc Yoshidal (|2004| ): fllennawi fc Sperge' 



l 2005h : lMaturi et al.1 ^200^ . 



rgil 



Let us first defined the weak lensing mass map which is 
the smoothed lensing convergence field (k): 



(10) 



/c(6i) = / d> K{4,-e)u{\4>\] 



where U is the filter function to be specified below. The same 
quantity is obtained from the shear data by 



(11) 



/C(6») = / d>7t(</):6»)Q(| 



where 7t(</) : 9) is the tangential component of the shear at 
position (j) relative to the point 9, and Q relates to U by 

Q{e)^ f de' 9'U{6') ~U{6). (12) 

^0 

We consider Q with a finite extent; in this case, one finds 

(13) 

where 9o is the outer boundary of the filter. Note that this 
is equivalent to set U a finite compensated filter; that is, 
j"" de 9 U{9) = and U{9) = for e > 

The basic idea of weak lensing cluster finding is to first 
construct a weak lensing mass map by applying eq. (|ll|l 
to shear data, then to search for high peaks in the map 
which are plausible candidates of massive clusters. The root- 
mean-square (rms) noise coming from intrinsic ellipticity of 
gal axies (whic h we call the galaxy shape noise) is evaluated 
by (|Schneider| [r996'). 



2 



2n„ 



d9 9 Q {9), 



(14) 



where <Je is the rms value of intrinsic ellipticityies of galax- 
ies and rig is the number density of galaxies. Throughout 
this paper, we take = 0.4 and rig = 30 arcmin"^, which 
resembles the shape noise expected for the Hyper Suprime- 
Cam survey. The signal-to-noise ratio (S/N) of the weak 
lensing map is defined by the ratio between the peak height 
and (Tshapc, 

^ = (15) 

^shapc 

In this paper we consider the following three filter func- 
tions. One is the truncated Gaussian (for U), 



Qcie) 



1 



1-1 + 



exp 



(16) 



for 9 < 9o and Qg = elsew here. We take 9g = 1' , which 
is the value also adopted in iHamana. Takada. fc Yoshidal 
(|2004h . and 9o = 15'. The others filters have the follow- 
ing functional form, consisting of the power-law with outer 
exponential-cutoff. 



^^(l + a6»/e»/){2+") 



exp 



291 



(17) 
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Figure 1. Three window functions considered in this paper are 
plotted. Top left: Absolute value of U is shown, with negative 
U being plotted by thin lines. The normalization is arbitrary. 
The convergence profile of the smoothl y truncated NF\y halo 
(Mvir = W^h-^MQ, zi = 0.3, Zs = 1) jOeuri fc Hamanall2011^ 
is also plotted as the dotted line. Bottom left: 9^Kjqp\fi/U, which 
is the contribution to the K from the logarithmic 8 interval, is 
plotted. Top right: The filter Q. The normalization is arbitrary. 
The tangential shear profile of the smoothly truncated NFW halo 
is also plotted. Bottom right: B^^tQ, which is the contribution to 
the K from the logarithmic 9 interval, is plotted. 



We consider two cases: (n, a) = (0,0.25) and (1,0.7) which 
we call PEXO and PEXl, res pectively. The form er mimics 
the filter function proposed bv lMaturi et al.l (|2005l l designed 
for maximizing the S/N of weak lensing peak by the NFW 
halo relative to noises coming from galaxy shape and cos- 
mi c structures, wh ereas the latter mimics the one proposed 
bv lHennawi fc s"pergel (2005.) which has a similar shape at 
the outer region but has less power on the inner region 
for suppressing the galaxy shape noise. We take Of = l', 
9e = 5' and 9o = 15' for both the cases, which are cho- 
sen so as to maximize the S/N for a cluster at z = 0.3 with 
Mvir = lO^*/i~^M0. Note that the filter scales do not need to 
be fi xed but in general can be varied to act as a matched filter 
(e.g., iHennawi k, Spergeill2005l : iMarian. Smith, fc BernsteinI 
I2OIOI I. In this paper we do not vary the filter scales because 
it is beyond the scope of this paper. The shapes of the filter 
functions are plotted in two top panels of Fig. [T] 

The weak lensing peak height for the NFW halo is com- 
puted by 



/Cnfw = 27r / d4> 4> K-MFw{4>)U{4>), 



(18) 



where knfw is the convergence profile from NFW halo for 
which we take the sm o othly truncated N FW profile (see 
lOguri fc Hamanal I2OIII : iBaltz et al] |2009| . for analytic ex- 
pressions). It is seen from bottom panels of Fig. [1] that the 
most of the contribution to /Cnfw comes from the mat- 
ter within the scale radius or from shear data on scales 
Os < e < Svir. We denote the S/N expected from NFW 
halo by !/nfw ~ /CNFw/o"ahape- The expected S/N is shown 
in Fig. [2] as a function of the concentration parameter. We 
find that, in the cases of lower A, the S/N is sensitive to 
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Figure 2. Top: The peak S/N expected for NFW halo as a func- 
tion of the concentration parameter. Different lines show different 
spherical overdensity mass shown in unit of Vd^'^h-'^ Mq. The lens 
and source redshift are set to be 0.3 and 1, respectively. Bottom: 
Corresponding virial mass is plotted. 



the concentration paramet er (the sa me argument but for the 
case of Mvir was made by iKing fc Mead (2011)). This can 
be explained by the fact that the lensing S/N mostly comes 
from the mass at the inner region of the halo as was shown in 
Fig. [1] The change in the concentration parameter results in 
the change in the mass of the inner region, and thus results 
in the change in the S/N . On the other hand, in the cases of 
higher A, the spherical overdensity mass effectively defines 
the mass of inner region; therefore, the change in the concen- 
tration parameter does not affect the S/N , though it alters 
the virial mass. Therefore we can argue that the peak S/N is 
not a good virial mass indicator but more tightly correlates 
with the inner mass such as M iqoq . The result app ears to be 
consistent with the finding of lOkabe et al.l l|2010l) . who ar- 
gued that weak lensing mass measurement errors are smaller 
for larger overdensities of A = 500 — 2000 than the case for 
the virial overdensity. 



3 NUMERICAL SIMULATIONS 

3.1 Gravitational lensing ray-tracing simulations 

We use a large set of ray-t racin g simulations that are de - 
tailed in ISato et al.l Hooi) and lOguri fc Hamanal l|201lh . 
It is based on the ray -tracing technique developed in 
iHamana fc Mellieil l|2001 In what follows we describe only 
aspects directly relevant to this study. 

The ray-tracing simulations are based on 200 real- 
izations of Af-body simulations with the box sizes of 

We 



240/i"^ Mpc (particle mass = 5.4 x lO^^/i'^Me 



^ The ray-tracing simulation codes are publicly available at 
|http : //th ■ nao ■ ac . j p/MEMBER/hamanatk/RAYTRIX/ Index ■ html [ 
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use the standard multiple lens plane algorithm to simulate 
gravitational lensing by intervening matter. In this study, 
we consider a single source redshift of Zs = 1, which is 
a typical mean source redshift of weak lensing analysis. 
Note that if redshift information on the source galaxies, 
e.g. from photometric redshift, is available, one may em- 
ploy the topographic technique which i mproves the capabil- 
ity of weak lensing cluster identification (jHennawi &: Spergell 
I2OO5I : [pTetrich fc Hartlap 2010 ). Instead of using 1 000 ray- 
tracing realizations generated by ISato et akl (|2009l ) , we re- 
generated 200 independent realizations so that the same halo 
does not appear in different realizations. 

Weak lensing mass maps, /C, are generated from lens- 
ing shear data by applying the operation eq. (Illf) on grid 
points of 2048 x 2048 with the grid spacing of 0.15 arcmin. 
We generated two types of mass maps, one is maps with- 
out the galaxy shape noise, and the other is maps with the 
galaxy shape noise. For galaxy shape noise, random ellip- 
ticities drawn from the truncated two-dimensional Gaussian 
(see below) are added to shear data: 

with CTe = 0.4 and assumed number density of source galax- 
ies of 30 arcmin"^. Peaks are identified as pixels that have 
higher values of K, than 8 surrounding pixels. Here it should 
be noted that in the case of the noise added maps, especially 
for the case of the PEXO filter, crowds of high peaks caused 
by the galaxy shape noise tends to form at massive cluster 
regions. To exclude duplicated detections in each cluster, we 
apply friends-of-fri ends (FOF) algorit hm with the linking 
length of 1 arcmin (|Maturi et al.|[2005h . and group together 
duplicated peaks as a single peak. This procedure generates 
a catalogue of lensing peaks for each map. 

3.2 Dark matter halo catalogue 

For each A'^-body output, we identify dark matter haloes 
using the standard FOF algorithm with the linking param- 
eter of fe = 0.2 and derive a halo mass Afpop for each halo. 
We define the halo position by the potential minimum of 
haloes where the gravitational potential is computed from 
only FOF member particles. In addition to Mpop, we com- 
pute the spherical overdensity mass with average overdensi- 
ties of A = 500 and 1000, which are denoted by M500 and 
A^iooo, respectively. 

To evaluate the halo shape, we compute the inertia ten- 
sor of the mass distribution: 

iVp 

= m-p'^ixi^k - Xi){xj^k - Xj), (20) 

k=l 

where i,j run from 1 to 3, A'j, is the number of FOF member 
particles of haloes and Xi denotes the cluster centre defined 
by the potential minimum. We use the trace of the inertia 
tensor to estimate the concentration of the halo mass distri- 
bution, 

7= i(/ll-h/22 +/33). (21) 



I I I I I . I I 'I . I ' I I. 




Figure 3. Peak values in the noise-free K map, u^, is compared 
with the NFW-corresponding peak S/N computed from Afpop 
{bottom left), M500 {iop left) and Miooo {top right). The mass 
maps generated from the Gaussian filter are used. Each point 
shows the relation for each halo, and the solid line shows the 
mean of them. Here haloes with Afpop > 1 X Mq at 

O.l < z < 0.4 are considered. 

This is compared with the expected value for the uniform 
overdens sphere: 

/th = ^A^virr^i,. (22) 
5 

To evaluate the triaxility of halo shape, we compute the 
eigenvectors from the inertia tensor and convert them to the 
axial ratios a/c and a/b {a ^ b ^ c). We also compute the 
angle between the line of sight and the major-axis direction 
of the halo which we denote as 6^ ■ 

Halo catalogues on the light cone are generated by 
stacking the simulation outputs in the same manner as in 
the ray-tracing experiments. In summary, each halo in the 
mock catalogues has data on the mass, redshift, the angular 
position on the weak lensing mass map and the parame- 
ters for the halo shape I, a/c and 6^. In addition, for each 
halo, we compute the NFW-corresponding peak S/N using 
eq. H18p adopting the mass-concentration relationship eq. 
for three mass estimates, Afpop, Afsoo and Afiooo, which 
we denote i^fof, z^soo and i^iooo respectively. Also we com- 
pute the virial radius defined by eq. and 7th defined by 
eq. (|22p . adopting Afpop as the virial mass. 



4 RESULT 

4.1 Halo-peak relationship 

We start with matching halo catalogues with weak lensing 
peaks. To do so, for each halo, the highest peak within the 
virial radius of the halo (or within 3 arcmin if the virial ra- 
dius is larger than 3 arcmin) is searched for and is considered 
as the corresponding peak. We denote the peak heights by 
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Figure 4. Relationship between viooo and Vk as a function of 
1^1000 ■ Tlie dots sliow tlie relation for each halo and the filled 
circles with error bars show the mean and the rms scatter, and 
the dashed curve show the value of the rms scatter. The K, maps 
generated from the Gaussian filter are used, and only the haloes 
with Miooo > 3 X lOl^?i-^Af0 at 0.1 < z < 0.4 arc considered. 
The top and bottom panels show results without and with the 
shape noise, respectively. 



First we examine the correlation between the Uk from 
the noise-free K map and the NFW-corresponding peaks 
evaluated from the three definitions of the halo mass, Mfof, 
M500 and Miooo- The aim of this analysis is to check the 
finding in Section 12.21 that the mass definition taking into 
account only the inner region (such as Miooo) better corre- 
lates with the peak S/N than masses including outer regions 
(such as Af2oo or Mvir). Fig. [3] compares Uk. with three NFW- 
corresponding peaks J^oF, v^oa and J/iooo, where haloes with 
MpoF > 1 X 10^''/i"^Mq at Q.l < z < 0.4 are considered. 
To quantify the correlation, we evaluate the mean and rms 
scatter of 5v = {vk — Vi)lvi among peaks with Vi > 4, 
and we find (mean, RMS)= (-0.11,0.25), (0.017,0.20) and 
(0.0098, 0.19) for i = FOF, 500 and 1000, respectively. Thus 
we find that i^iooo best correlates with Vk., as expected. 
Based on this result, in what follows, we take Miooo as 
the mass indicator of the haloes. Note that for the NFW 
halo with Cvir = 4, Miooo/Mvir — 0.35, and Miooo /Afvir 
is larger (smaller) for a halo with higher (lower) concen- 
tration. We also find that vyov-v^ relationships not only 
have larger scatters but also show some offset in the mean 
relationship. One of the reasons may be the mismatch be- 
tween Mfof computed in A^- body simulations and M yir de- 
fined by overdensity (see also lOguri fc Ifamanallioill ). pre- 
sumably caused by the presence of outer substructures that 
contribute to the FOF mass. In the analysis above, the mass 
maps generated from the Gaussian filter are used though the 
results are found not to be affected by choice of the filter. 

Next we examine the influence of the following three 



effects on the peak SjN; (1) the galaxy shape noise, i.e., 
the noise due to the intrinsic shapes of source galaxies used 
for weak lensing analysis, (2) the projection of structures 
along the line of sight of clusters, which we call the cosmic 
noise, and (3) the diversity of dark matter distributions of 
haloes, which we call the halo shape effect. To do so, the mass 
maps generated from the Gaussian filter are used though the 
results are found not to be strongly dependent on the choice 
of the filter. With our choice of parameters for the galaxy 
shape noise (e.g., ae — 0.4 and rig — 30 arcmin^^), the rms 
of the galaxy shape noise is cTshapo = 0.02, which should act 
as the additive noise in The a nalytic wa y to e stimate 
the cosmic noise was developed by llfoekstral (|200lh under 
the assumption that the halo and the large-scale structures 
are uncorrelated, and for the cosmological model adopted 
in this paper and the Gaussian window function, its rms 
is estimated to be acosmic = 0.01, which should again act 
as the additive noise. The halo shape effect should, in a 
crude approximation, act as the multiplicative noise in the 
peak S/m 

We estimate its amplitude using the ray-tracing 
simulation data. We evaluate scatters in j/„ — i/iooo relations 
as a function of z^iooo , and results are plotted in Fig. |4l 

The result for the noise-free IC map case shown in the 
top panel of Fig. |4] indicates that the rms increase with 
i^iooo- Since only the cosmic noise and halo shape effect are 
playing in this shape noise free case, and the cosmic noise 
should be independent of the S/N, the dependence of the 
rms on !/iooo is most likely originated from the halo shape 
effect. We estimate from the rms measured from — i^iooo 
relationship and the expectation value of acosmic = 0.01 
(thus acosmic/ashapo = 0.5) that the rms of the halo shape 
effect scales approximately ahaio ~ (0.1 — 0.2) x JCnfw 
(i.e., ahaio/ashapc ~ (0.1 - 0.2) x u). Thus for very mas- 
sive haloes the halo shape effects can be larger than the 
other noises (the q ualitatively same argument was made by 
iTan fc FanI (|2005h '). However, it should be noted that for 
a very shallow survey the above argument may not be the 
case as ashapc scales with the source galaxy number density 
as ashapc x nj^'^. In the noise added case plotted in the 
bottom panel of Fig. [J] the mean value is slightly greater 
than zero, quantitatively {v^ — J'looo) ~ 0.25 thus ~ 5 per 
cent for the v = 5 peak. This is because we consider the 
peaks that are affec ted by the shape noise in an a s ymme tric 
way as discussed in iHamana. Takada. fc Yoshidal (|2004h . 

We now discuss details of the halo shape effect. In Fig.O 
the fractional difference between I'looo and Vk measured from 
IC map (with the galaxy shape noise added) is shown as a 
function of the halo concentration I/Ith (top-left panel), the 
axis ratio a/c (bottom-left panel) and the halo orientation 
with respect to the line-of-sight direction |cos(^2)| (right 
panels). Dots indicate values for each halo, where we con- 
sider haloes with Miooo > 3 x W^^H'^Mq at 0.1 < z < 0.4 
with the mean S/N of this halo sample being (t/iooo) = 3.7. 



^ Approximately, deviations of halo shapes from a fiducial spher- 
ical model can be written as p ~ (l + '5)Psph where 5 describes de- 
viations from the fiducial spherical mass distribution denoted by 
by Paph. The lensing convergence, k, is the line-of-sight projection 
of the halo mass distribution weighted by the lensing efficiency, 
and thus can be approximated as k ~ (l-l-5)Kaph, indicating that 
the effect of halo shapes is multiplicative. 
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Figure 5. Relationship between J^iooo and as a function of 
I/Ith {top 'fi/Oi '^/'^ {bottom left) and |cos(S2)| {right panels). 
Here lialoes witti Miooo > 3 X IO^^/i'^Mq at 0.1 < z < 0.4 are 
considered. The dots show the relationship for each halo and filled 
circles show mean, vertical error bars show the rms scatter, and 
the horizontal error bars show the range where sample is taken. 
The K maps with the shape noise generated from the Gaussian 
filter are used except for the bottom-right panel where the noise- 
free K. maps are used. The green curve in the bottom-right panel 
shows the theoretical prediction based on the triaxial halo model. 



Filled circles show average values. We find that there ex- 
ists (1) a clear correlation between the halo orientation and 
the peak height deviations from the NFW model prediction, 
whereas (2) no correlation for the halo shape parameters 
{I/Ith and a/c). The results indicate that the intrinsic halo 
shape (concentration and axis ratio) does not cause a sys- 
tematic bias in the peak heights as long as one employs 
an appropriate definition of the halo mass such as A'/iooo, 
but just contribute to the scatter. However the halo orienta- 
tion does cause the systematic bias because the line-of-sight 
projected mass at the inn er region depends strongly on th e 
halo orientation (see also lOguri et al.[|2005l : lGavazzill2005l ). 
In the bottom-right panel of Fig. [SI the theoretical predic- 
tion based on the triajci al halo model (Jing & Suto 2002; 
lOguri fc BlandfordI |2009| ) is also shown (see Section l2.ip . 
We find that the triaxial model nicely reproduces the orien- 
tation dependence of the peak heights found in ray-tracing 
simulations, except for the small offset which is originated 
from an approxim ation involved in the triaxial haloes (see 
lOguri et al.[|2005l ). 
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Figure 6. Probability distribution functions of the halo orien- 
tation with respect to the line-of-sight direction, |cos(6z)[, for 
weak lensing selected clusters. We adopt the peak height thresh- 
old value of u = 4. The dashed histogram is for the NFW- 
corresponding peaks, i.e., peak heights evaluated by the spherical 
NFW model via spherical overdensity masses, the long-dashed 
histogram for the peak heights measured from the noise-free /C 
map, and the solid histogram for the peak heights measured from 
the noise added K map. The dotted curve shows the analytic 
model expectation based on the triaxial halo model. 
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Figure 7. The halo orientation bias as a function of the peak 
S/N threshold. We quantify the bias by the ratio of the numbers 
of haloes with the major axis aligned to the line-of-sight direction 
{\cos8z\ > 0.5) to those with anti-aligned (|cos9z| < 0.5). The 
dashed histogram is for the peak heights measured from the noise- 
free /C map, and the solid histogram for the peak heights measured 
from the noise added K, map. The dotted histogram shows the 
analytic model prediction based on the triaxial halo model. 



4.2 Selection bias in weak lensing selected clusters 

Given the strong dependence of the peak heights on the halo 
orientation, we shall now investigate its impact on weak 
lensing selected cluster catalogues. In Fig. [51 we show the 
probability distribution function (PDF) of the halo orienta- 
tion with respect to the line-of-sight direction, | cos{9z)\, for 
haloes with peak heights above the threshold value u = 4. 
Here we include only haloes at 0.1 < z < 0.4 for sim- 
plicity. Different histograms are for different causes; NFW- 



corresponding peaks (dashed), peak heights measured from 
the noise- free /C map (long-dashed), and the noise added 
map (solid). In addition, the same PDF computed with 
analytic triaxial halo model (see Section 12.11) is also plot- 
ted by the dotted line, which should be compared with the 
noise-free case because the galaxy shape noise is not added 
in the triaxial model calculation. The PDF for the NFW- 
corresponding peaks (i^iooo) is flat. This is exactly expected 
because the NFW-corresponding peak heights are computed 
from Miooo assuming the spherical NFW profile, and hence 
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they should not depend on the orientations of haloes. On 
the other hand, the PDFs for the measured peaks are signif- 
icantly skewed such that the major axes of haloes are prefer- 
entially aligned with the line-of-sight direction. To quantify 
the bias, we introduce an estimator defined by the ratio be- 
tween the numbers of haloes with the major axis aligned to 
the line-of-sight direction (IcosS^j > 0.5) and anti-aligned 
(|cos6'^| < 0.5), i.e., = n{\cos9z\ > 0.5)/n{\cos0^\ < 
0.5). For the cases shown in Fig.[6l the bias are = 1.7 and 
1.4 for the noise- free and noise added cases, respectively. For 
the case of the NFW-corresponding peaks we have = 1.07, 
suggesting that the error in the above estimation is less than 
10 per cent. 

In addition, we examine the bias as a function of the 
peak heights using the estimator, hzii^n) = ti{vk.', Icos^^I > 
0.5)/n(//K; |cos6'z| < 0.5). The results are shown in Fig. [T] 
While the statistic is not very good especially for higher S/N 
bins because of a limited number of peaks, one can clearly 
see a trend that the orientation bias is larger for higher S/N . 
This is presumably because the halo shape effect is multi- 
plicative to the S/N . Although the results are not accurate 
enough to make precise predictions, we find that the bias is 
^ 1.5 for high S/N peaks of > 5. and ~ 1.2 for ~ 4. 
Therefore the orientation bias is one of the most significant 
selection bias in weak lensing selected cluster catalogues. 
Its impact on cluster related sciences should be taken into 
consideration. The results shown in Figs. [3 [6l and [7] indi- 
cates that the analytic triax;ial model predictions are in very 
good agreement with the trends found in ray-tracing simu- 
lations. This confirms that the orientation bias comes from 
the intrinsic triaxial shapes of haloes, rather than the effect 
of surrounding matter around clusters, because such corre- 
lated matter is not included in our analytic triaxial model 
calculations. 

We note that the strong halo orientation bias was 
also found in previ ous strong lensing studies as well. 
iHennawi et al.l (|2007l l used a large set of ray-tracing sim- 
ulations to show that major axes of clusters producing giant 
arcs are preferentially aligned with the line-of-sight direc- 
tion, with the level of the bias si milar to what found in this 
paper. [Oguri &: BlandfordI l|2009l ) reported that clusters hav- 
ing larger Einstein radii have larger orientation bias, based 
on the triaixl halo model. Si milar result for st r ong le nsing 
clusters was obtained also bv'Meneg hetti et al.l (|2010l ). Our 
result indicates that similar orientation bias exists for weak 
lensing selected clusters as well. 

4.3 Completeness 

To see the impact of the halo orientation bias, we exam- 
ine the completeness which we define by the fraction of 
haloes identified as high peaks above a given threshold 
{vmin) in the /C map relative to all the haloes, namely, 
n{M, z; Vk > i'min)/n{M, z). The choice of the threshold is 
arbitrary, usually it is chosen considering the trade-off be- 
tween the sample size and the fraction of the false positives 
(see Section [4.4p . Here we adopt fmin = 4. The complete- 
ness for different filters are shown in left panels of Fig. [S] We 
find that the PEXO filter is the best in terms of the com- 
pleteness, although the difference between filters is small, 
10 per cent at most. Right panels of Fig. [8] compares the 
completeness for haloes with IcosS^] > 0.5 (red histogram) 




Figure 8. Completeness as a function of the halo mass. From 
top to bottom panels, we change the redshift range from lower 
to higher mean redshifts. Left panels: Different histograms show 
results for different filters. Blue, red, and black histograms are for 
the PEXO, PSXl and Gaussian filters, respectively. Right panels: 
The red histogram is for the the haloes with |cos(92)| > 0.5, i.e., 
haloes with major axes aligned with the line-of-sight direction, 
whereas the blue histogram is for anti-aligned haloes. 



Table 1. Summary of the number density (per 1 deg^) of weak 
lensing mass peaks (in noise added K maps) above thresholds. 



filter 


Npc^kii^K > 4) 


Npc^kii^K > 5) 


Npe^U^. > 6) 


PEXO 


3.3 


0.70 


0.22 


PEXl 


2.7 


0.59 


0.19 


Gauss 


2.2 


0.53 


0.19 



with the others (blue), where the IC maps with PEXO fil- 
ter are used. We find that the effect of the halo orientation 
bias is visible for the intermediate mass haloes for which 
z/NFW ~ J^min. The reason is that for the very massive/small 
haloes, the S/N differs very much from i^min, meaning that 
the spread of v produced by the halo shape effect does not 
affect the detectability of those haloes. Thus the complete- 
ness of haloes with z/nfw ~ i^min are mostly affected by the 
halo orientation bias. The difference of the amplitude for the 
completeness is about 20 per cent for those haloes. 

4.4 Peak counts and purity 

So far we discussed connections between the weak lensing 
properties of clusters and dark matter halo properties. Here, 
we study the subject from a different viewpoint, i.e., we 
examine properties of the weak lensing peak sample in terms 
of peak counts and purity, paying special attention to the 
comparison among three filters. 

First we examine peak counts and compare with the 
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Figure 9. Number counts of weak lensing peaks for three fil- 
ters from top to bottom panels. The red and blue histograms 
are measurements from ray-tracing simulations with and with- 
out galaxy shape noise, respectively. The green histogram is the 
number count of peaks from the pure noise map. The long-dashed 
and dashed curves are analytical predictions using the spherical 
NFW profile, with and without shape noise taken into account, 
respectively. 



analytic model. In Fig. peak counts measured from three 
K. maps (the galaxy noise free, noise added, and the pure 
noise cases) are presented for three filters. The number 
densities of peaks (in noise added JC maps) above the 
given thresholds are summarized in Table [T] We find that 
the PEXO filter yields the largest peak counts, specifically 
about 20 per cent larger than the PEXl filter, and even 
larger difference from the Gaussian filter. This is a nat- 
ural consequence that the PEXO filt e r m imics the opti- 
mal filter developed by iMaturi et al. I (12005 ) which is de- 
sign ed for maximizing the S/N of JC from the NFW halo 
(see IPace et"al] (|2007| 'I for the test of the capability of the 
optimal filter against numerical simulations). To compute 
the analytic prediction of the peak cou nts, we follows the 
so-cal led halo mode l first developed by Kruse fc Schneider] 
||2000D (see also iBartelniann. King fc Schneider! l200ll : 
iHamana. Takada. fc Yoshidall2004^ ■ in which it is assumed 
that high peaks are dominated by lensing signals from sin- 
gle massive haloes. Here we adopt the spherical NFW den- 
sity profile for the analytic calculation of the number count, 
because we include the effect of the halo triaxiality via 
the halo shape noise as described below. To take into ac- 
count the effect of no ise, we employ the approxim a te ap - 
proach developed by iHamana. Takada. fc Yoshidal (|2004l ') 
with some modification. Specifically, we make the follow- 
ing three assumptions. (1) Very high peaks are neither re- 
moved nor generated by the noise but their peak heights 
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Figure 10. The purity for the three filters are plotted by lines in 
different colors. The dashed histogram is for the noise-free case. 
Top panel is for the case where peaks are matched with haloes 
with Miooo > 3 X 10^'^h~^ Mq, whereas bottom panel is for 
A^lOOO > 1-5 X 10^^h~^MQ. The horizontal dotted line shows the 
probability that a peak is matched with a randomly distributed 
halo by chance. 



are altered by the noise. (2) The scatter in peak heights 
with respect to the corresponding NFW peak height fol- 
lows the Gaussian distribution with the standard dev iation 
of o-pcak- While in IHaman a. Taka da. fc Yoshidal (|2004h only 
the galaxy shape noise was taken into accounted, here we 
include both the cosmic noise and halo shape noise, in ad- 
dition to the galaxy shape noise. (3) Within a small range 
of the peak height (ly), the peak counts can be approxi- 
mated by the exponential form npcak(i') = exp(p!/), with 
a constant exponential index of p. Under these assumptions 
the peak counts in the presence of the noises are given by 
IHamana. Takada. fc Yoshidall2004l ). 



~ exp 



fV 



ripaakil^), 



(23) 



where / — cTpoak/cnoisc . To estimate /, we assume the pro- 
portion of (Tshapc : o-cosmic : o-haio = 2:1:1. For the case of 
the Gaussian filter, this corresponds to (Thaio = 0.01 which is 
a reasonable approximation for peaks with ~ 5 (see Sec- 
tion l4.Hl . We take the quadrature sum of these three compo- 
nents to obtain / ~ 1.2. The peak counts computed in this 
method are presented in Fig. [51 We find that for high S/N of 
I' > 5 the improved analytic prediction agrees well with the 
measurements, whereas for lower S/N the measured counts 
are slightly larger than the prediction. A possible origin of 
the excess is the false signals which we examine below. 

Finally, we examine the purity which we here define by 
the fraction of peaks which are associated with true haloes 
among all the peaks. However, in practice, the purity is not 
uniquely defined quantity because of the following two ambi- 
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guities. One is the allowed separation between peak position 
and halo position. The other is the minimum mass of haloes 
allowed to be associated with peaks. We adopt the minimum 
halo mass of Miooo = 3 x lO^^/i-^M© or 1.5 x lO^ft-^M© 
based on the detectability study shown in Fig. [S] Consid- 
ering the virial radii of those haloes, we fix the maximum 
separation at 3 arcmin. The results are presented in Fig. 1101 
We also plot the probability that a peak matched with a 
randomly distributed halo (with the same number density 
as the true halo catalogue) by chance by the horizontal dot- 
ted line, which gives an estimate of the chance coincidence 
between peaks and unrelated haloes. We find that the proba- 
bility of the chance matching is not significant for our choice 
of the minimum masses and the allowed separation. The pu- 
rity is not dependent on the choice of the filter, and is very 
high (i.e., the false positive rate is very small) for peaks with 
u > However the purity drops rapidly at lower S/N , and 
it becomes about 50 per cent at ~ 4. We argue that this 
accounts for the excess in the peak counts over the theoret- 
ical prediction found in Fig. [9] In the same figure, we also 
plot the purity measured from the galaxy shape noise free 
case. The false positives in this case may arise from chance 
projections such as the line-of-sight projection of small mul- 
tiple haloes or filamentary structure. As the purity for the 
noise-free case is found to be greater than 90 per cent for 
V > A, we conclude that the contamination of such chance 
projection is not significant. 



5 COMPARISON WITH OBSERVATIONS 

It is worth checking the results of our mock simulation data 
against real observational results. In Fig. 1111 we compare the 
peak counts measured from real weak lensing K maps gener- 
ated from 4 deg^ of Subaru/Suprime-Cam data with results 
from the mock simulations. A description of the Suprime- 
Cam data and data analysis is given in Appendix [X] In 
short, the Suprime-Cam data consist of four fields with fairly 
uniform data quality (with respect to the depth and see- 
ing condition) with the number density of galaxies used for 
weak lensing analysis is Ug ~ 27 arcmin"^ and the rms 
ellipticity of (Tshape = 0.4 which are in a good agreement 
with the values adopted in the mock simulations. The mean 
source redshift is not known as there is no large spectro- 
scopic/photometric redshift galaxy catalogue reaching the 
depth of the data {j! ~ 25.5 AB mag), but the value of 
(zs) = 1 assumed in mock simulations is reasonable (see also 
lOguri et al]|2012f ). We generated K. maps with the three fil- 
ters under consideration, and we searched for peaks in the 
/C maps. The peaks located within 1 arcmin from the field 
boundary were discarded as the regions are likely affected 
by the partial lack of data. The total area used for peak 
finding is 4 deg^ . We detect 14, 9 and 7 peaks with v > A for 
the PEXO, PEXl and Gaussian filter, respectively. To mea- 
sure the peak counts from mock simulations under a similar 
survey condition, we extracted a contiguous 4 deg^ region 
from each mock weak lensing realization, and measured peak 
counts. We computed the average counts and the range en- 
closing 68 per cent of 200 realizations and show in Fig. 1111 
We find that (i) peak counts from real data are in a reason- 
able agreement with the expectation from mock data, and 
(ii) its dependence on the filter is similar to that expected 




Figure 11. Red histogram shows the number count of peaks 
measured from real weak lensing K maps generated from 4 deg^ 
of Suprime-Cam data (see Appendix lAl for details). Results from 
mock simulation of 4 deg^ survey are shown by black histogram 
(the average among 200 realizations) with shading (the range en- 
closing 68 per cent). Note that on the highest bin, peak counts 
are zero in most of the realizations and non-zero counts are mea- 
sured in a small fraction of realizations, therefore the average is 
non-zero but the 68 per cent enclosing region is not defined. 

from the mock data. Of course this comparison tests only a 
limited aspect of the mock simulation, yet the agreements 
can be regarded as a piece of evidence that our numerical 
simulations based on the standard cosmological model and 
realistic noise parameters produce reasonably realistic mock 
data. 



6 SUMMARY AND DISCUSSIONS 

We have investigated scatter and bias in weak lensing se- 
lected clusters, employing both the analytic model descrip- 
tion of dark matter haloes and the numerical mock data 
of weak lensing cluster surveys generated with gravitational 
lensing ray-tracing through a large set of A'^-body simula- 
tions. We have paid special attention to the effects of di- 
versity in the dark matter distribution of haloes. Our major 
findings are summarized as follows. 

(1) We have examined the relationship between the peaks 
measured from the noise-free K, map and the expected peak 
values computed assuming the spherical NEW profile with 
three cluster mass definitions, Mfof, M500 and Miooo- We 
have found that the expected peak value computed with 
Miooo best correlates with the measured peak heights. This 
confirms the finding in Section 12.21 that the lensing S/N 
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computed with the spherical overdensity mass with higher 
mean interior density is less affected by variation in halo 
concentrations than those with lower interior density. An 
implication to observational studies is that the peak S/N is 
not a good virial mass indicator but is more tightly related 
to the inner mass such as Miooo- 

(2) We have examined the influence of the diversity of 
the halo shape on the peak S/N values by comparing the 
peaks measured from the noise-free K. map with the NFW- 
corresponding peak heights computed with Miooo using the 
spherical NFW profile. We have found that the rms of the 
halo shape effect scales approximately (Jhaio ~ (0.1 — 0.2) x 
K.NFW (i.e., ahaio/o-shapc ~ (0.1 - 0.2) X v). We therefore 
conclude that the scatter caused by the halo shape effect can 
be larger than the ones from other noises for very massive 
haloes. 

(3) We have found a clear correlation between the halo 
orientation and the peak heights, such that haloes whose ma- 
jor axes are aligned with the line-of-sight direction tend to 
generate a higher peak than the expectation of the spherical 
NFW model. Using both numerical and analytic approaches, 
we have examined the systematic bias caused by the orien- 
tation effect. We have evaluated the bias using the ratio 
between the numbers of haloes (identified with peak S/N 
above certain thresholds) with the major axis aligned with 
the line-of-sight direction (IcosS^I > 0.5) and those anti- 
aligned. We found that the bias is hz — 1.4 for the threshold 
of Vmin = 4, Indicating that orientation bias is quite strong. 
In addition, we have found that the bias is larger for higher 
S/N peaks, because the halo shape effect is multiplicative to 
S/N. Thus the orientation bias is a non-negligible selection 
bias in weak lensing selected cluster catalogues. We have 
also examined the effect of the halo orientation bias on the 
completeness and found that the haloes with vnfw ~ J^min 
are most affected by the halo orientation bias, and its ampli- 
tude for the completeness is found to be about 20 per cent 
for those haloes. We have shown that the analytic triaxial 
halo model explains the orientation bias found in ray-tracing 
simulations very well, which confirms that the halo triaxial- 
ity is the main cause of the orientation bias. 

(4) We have compared the capability of three filters with 
respect to the the completeness, peak counts and purity. We 
have found that the PEXO works best, which gives about 10 
per cent better completeness for haloes with f ^ Vmin than 
the other two filters and about > 20 per cent larger peak 
counts with the similar level of purity. 

(5) We have developed a prescription to an- 
alytically compute the number count of weak 
lensing peaks which improve earlier works 
dKruse fc Schneider] I2OO0I: iBartelrnann. King fc Schneider! 



I2OOII : iHamana. Takada. fc Yoshidal l2004h . Specificallv. we 
have employed an approximate way to include the scatters 
in peak heights caused by all the three sources addressed in 
this paper; the galaxy shape noise, the cosmic noise, and 
the halo shape noise. We have tested the improved model 
against the mock ray-tracing simulation data, and found 
that the analytic predictions agree well with the simulation 
results for a high S/N of J/ > 5. 

(6) We have compared the number counts of peaks from 
the numerical mock data with those from 4 deg'^ of real 
Subaru/Suprime-cam data. Although the comparison is lim- 
ited by small number statistics, we have found that both the 



number counts and their smoothing filter dependence are in 
reasonable agreement between the mock data and the ob- 
servation. 
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Table Al. Summary of data sets. See Figs. IA1IA4I for the coor- 
dinate of eacli field 



field name 
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rp b 

^ exp 


seeing'^ 




SXDS 


0.99 (0.77) 


6000 


0.52 


27.6 


COSMOS 


1.85 (1.32) 


2400 


0.53 


26.5 


Lockman-hole 


0.74 (0.56) 


3600 


0.51 


27.3 


ELAIS Nl 


1.99 (1.39) 


3600 


0.57 


26.6 



" Area after masking regions affected by bright stars in unit of 
deg^. The numbers in parentheses are the area after cutting the 
edge regions within 1 arcmin from the boundary. 
^ Exposure time for each pointing in the unit of seconds. 

Median of seeing full width at half-maximum in the unit of 
arcsec. 

Number density of galaxies used for the weak lensing analysis 
in the unit of arcmin"^. 
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APPENDIX A: WEAK LENSING ANALYSIS OF 
SUPRIME-CAM DATA 

We collected i'-band data taken with the Subaru/Suprime- 
Cam (Miyazaki et al. 2002) from the data archive SMOKj^ 
under the following three conditions: data are contiguous 
with at least four pointings, the exposure time for each 
pointing is longer than 1800 sec, and the seeing full width 
at half-maximum (FWHM) is better than 0.65 arcsec. Four 
data sets meet these requirements and are summarized in 
Table lAll 

E ach GCD data was reduced using the SDFre^ soft- 
ware (jYagi et all l2002l : lOuchi et ail l2004l l. Note that we 
conservatively use the data only within 15 arcmin ra- 
dius from the field centre of Suprime-Cam, because at 
the outside of that, the point spread function (PSF) be- 
comes elongated significantly which may make the correc- 
tion for the PSF inaccurate. Then mosaic s tacking was 
done with SCAMP (|Bertinll2006l ) and SWarjE^ (iBertin et all 
I2OO2I ). Object detections were performed with SExtractor 



* 'http : //smoka.nao . ac ■ jp/| 

^ In the process of the correction for both the field distor- 
tion and differential atmospheric dispersion, the bi-cubic resam- 

! iling scheme was inipleme nted to suppress the aliasing effect 
Hamana fc Mivazakill2008l) . 

SWarp was modified so that it can treat the bad pixel flag 
from the SDFred software properly. 
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Table A2. Summary of high peaks. 



fiolrl TTl 


R.A. 

[deg] 


Dec 

laegj 




(Peak-S/Af) 
PEXl GAUSS 


SXDS-1 


34.0185 


-5.0755 


4.1 


o.o 


O.O 


COSMOS-1 


150.5067 


2.7657 


5.2 


o.z 


0.1 


COSMOS-2 


150.2063 


1.8283 


4.3 


4.4 


4.4 


COSMOS-3 


150.1838 


1.6583 


4.0 


3.9 


3.7 


COSMOS-4 


150.0987 


2.7158 


4.0 


3.7 


2.9 


COSMOS-5 


149.8486 


2.2883 


5.1 


4.9 


4.3 


ELAIS-1 


243.3380 


55.4445 


4.0 


4.0 


3.8 


ELAIS-2 


242.8025 


55.5614 


6.5 


6.5 


6.6 


ELAIS-3 


242.6855 


55.4296 


4.2 


4.2 


3.6 


ELAIS-4 


241.2850 


55.6263 


4.2 


4.2 


4.2 


ELAIS-5 


241.2278 


54.5485 


4.2 


3.9 


3.4 


ELAIS-6 


241.2034 


54.4808 


4.1 


3.9 


3.4 


ELAIS-7 


241.1764 


54.5356 


5.3 


5.2 


5.2 


ELAIS-8 


240.8272 


55.1126 


4.3 


4.3 


4.3 



jBertin fc Arnoutd Il996l') and hUndpeaks of IMC AT soft- 
ware ( Kaiser. Squires fc BroadhurstI ll9S?5l '). and two cata- 
logues were merged by matching positions with a tolerance 
of 1 arcsec. 

For weak lensing measur ements, we follow the so-called 
KSB m ethod described in iKaiser. Squires fc Broadhurst 



199S 



1995 ). iLuppino fc Kaiser! (1997), and Hoekstra et al 



Stars are selected in a standard way by identify- 
ing the appropriate branch in the magnitude half-light ra- 
dius (rji) plane, along with the detection significance cut 
5'/A'^ > 10. The number density of stars is found to be 
~ 1 arcmin"^ for the four fields. We only use galaxies met 
the following three conditions, (i) the detection significance 
of S/N > 3 and nu > 10 where nu is an estimate of the peak 
significance given by hfindpeaks, (ii) is larger than the 
stellar branch, and (iii) the AB magnitude is in the range of 
22 <i' < 25.5 (where MAG_AUTO given by the SExtractor 
is used for the magnitude). The number density of resulting 
galaxy catalogue is quite uniform not only over each field 
but also among four fields (see Table [Alj). We measure the 
shapes of objects by getshapes of IMC AT, and correct for 
the PSF by the KSB method. The rms of the galaxy ellip- 
ticities after the PSF correction is found to be 0.4 for all the 
four fields. 

Weak lensing K, are computed from galaxy ellipticity 
(after the PSF correction) data with three filters under 
consideration on regular grid points with a grid spacing of 
0.15 arcmin. We then identify high peaks from the K maps. 
Peaks located within 1 arcmin from the field boundary were 
discarded because such regions are likely affected by the par- 
tial lack of data. The total area used for the peak finding 
is 4 deg^. We detect 14, 9 and 7 peaks with v > A for the 
PEXO, PEXl and Gaussian filter, respectively. High peaks 
with u > A are summarized in Table IA2I and peak counts 
are plotted in Fig. Illl in which a good agreement with mock 
simulation data is found. In Figs. IA1IIA4I K, maps gener- 
ated with the PEXO filter are shown in which high peaks 
are marked with circles. We also mark known clusters taken 
from a compilation by NASA/IPAC Extragalactic Database 
by the plus symbols. It is seen that many but not 




I 



34.5 
R.A. [degree] 

Figure Al. K maps for the SXDS field, generated with the PEXO 
filter are plotted by blue scale. Gray regions are either the region 
outside the analyzed area (within 15 arcmin from the field centre 
of Suprime-Cam) or masked regions where the data are affected 
by bright stars. High peaks with v > i are marked with red circle 
along with ID given in Table IA2I The Plus symbols show posi- 
tions of known clusters taken from a compilation by NASA /IPAC 
Extragalactic Database (NED). 
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Figure A2. Same as Fig.|XT]but for the COSMOS field. 

all high peaks are associated with known clusters. It should 
be noted that the known cluster sample plotted in the figures 
is just a compilation of many individual cluster catalogues 
identified by many different techniques, and thus is not ho- 
mogeneous even over a single field. 



|http : //ned . ipac ■ caltech . eduT] 
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Figure A3. Same as Fig. lAll but for ttie Lockman-liole iield. 
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Figure A4. Same as Fig.|AT]but for the ELAIS Nl field. 



